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SUMMARY 

Theoretical values of performance parameters for liquid ammonia and 
liquid fluorine as a rocket propellant were calculated on the assumption 
of equilibrium composition during the expansion process for a wide range 
of fuel-oxidant and expansion ratios. The parameters included were 
specific impulse, combustion-chamber temperature, nozzle-exit tenpera- 
ture, equilibrium composition, mean molecular weight, characteristic 
velocity, coefficient of thrust, ratio of nozzle-exit area to throat 
area, specific heat at constant pressure, coefficient of viscosity, and 
coefficient of thermal conductivity. 

The max imum value of specific impulse was 311.5 pound-seconds per 
pound for a chamber pressure of 300 pounds per square inch absolute 
(20.41 atm) and an exit pressure of 1 atmosphere. 


INTRODUCTION 

Liquid ammonia and liquid fluorine are of interest as a rocket pro- 
pellant because of high performance. Extensive data“exist in the litera- 
ture on their availability and cost, and on their physical, chemical, and 
handling properties. 

The perfor man ce of liquid ammonia and liquid fluorine has been 
reported in the literature by a number of organizations, such as Jet 
Propulsion Lab., C.I.T.; Reaction Motors, Inc.; Battelle Memorial 
Institute; and the MCA. Additional performance calculations for this 
propellant were made at the NACA Lewis laboratory as part of a series of 
calculations on propellants containing the chemical elements hydrogen, 
fluorine, and nitrogen (ref. l) to provide a comparison with the per- 
formance of other propellants based on the same thermodynamic data and 
computed to the same degree of accuracy, and to provide several parame- 
ters not previously published. 

Data were calculated on the basis of equilibrium composition during 
expansion and cover a wide range of fuel-o’xidant and expansion ratios. 
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The performance parameters included are specific impulse,' combustion- 
chamber temperature , nozzle -exit temperature, equilibrium composition, 
mean molecular weight^ characteristic velocity, coefficient of thrust, 
ratio of nozzle-exit area to throat area, specific heat at constant pres 
sure, coefficient of viscosity, and coefficient of thermal conductivity. 

So that data based on the assumptions of equilibrium and frozen com 
position during the expansion process could be compared, several addi- 
tional calculations were made assuming frozen composition. 


SYMBOLS 

The following symbols are used in this report: 

A number of equivalent formulas (a function of pressure and molecular 
weight; see ref. 5) 

a local velocity of sound, ft/ sec 

Cp ' coefficient of thrust 

C° molar specific heat at constant pressure, cal/ (mole ) (°K) 

Cp specific heat at constant pressure, cal/(gm)(°K) 

c specific heat at constant volume, cal/(gm)(°K) 

c* characteristic velocity, ft/sec 

(b log A \ 

A log T ) 

s 

(b log p ± \ 

\d log T ) s 

g acceleration due to gravity, 32.174 ft/sec^ 

Bp. sum of sensible enthalpy and chemical energy, cal/mole 

h sum of sensible enthalpy and chemical energy per unit weight, 

i 


nM 


cal/ g 
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specific impulse , lb-sec/lb 

coefficient of thermal conductivity, cal/(sec) (cm) (°K) 
molecular weight, g/mole 
number of moles 
pressure 
partial pressure 

universal gas constant (consistent units) 

equivalence ratio, ratio of number of fluorine atoms to hydro- 
gen atoms 

nozzle area, sq ft 

temperature, °K 

rate of flow, lb/sec 

d log a\ 
d log T/p 

3 log nA 
b log T/p 

d log P\ 
d log p j s 

p. coefficient of viscosity, g/ (cm) (sec) = poise 

p density, g/cu cm 

Subscripts : 

c combustion chamber 

e nozzle exit 

frozen composition assumed frozen 
i product of combustion 


max 


maximum 
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P constant pressure 

s constant entropy 

t nozzle throat 

x any point in nozzle 


CALCULATION OF PERFORMANCE DATA 

Calculations of the performance data were made with a Bell computer 
and an IBM Card-Programmed Electronic Calculator as described in refer- 
ence 1. The assumptions, thermodynamic data, and transport properties 
used for the calculations are the same as those of reference 1. 

The products of combustion were assumed to be ideal gases and 
included the following substances: hydrogen fluoride HF, hydrogen Hg, 

nitrogen Ng, fluorine Fg, atomic fluorine F, atomic hydrogen H, and 
atomic nitrogen N. The dissociation energy of Fg was taken to be 

35.6 kilocalories per mole (ref. 2). Physical and thermochemical proper- 
ties of the propellants were taken from references 2 to 5 and are given 
in table I. 

Procedure for combustion conditions . - For each of ten equivalence 
ratios, combustion temperature , equilibrium composition, enthalpy, mean 
molecular weight, derivative of the logarithm of pressure with respect to 
the logarithm of density at constant entropy y g , specific heat at con- 
stant pressure, coefficient of viscosity, coefficient of ther mal conduc- 
tivity, and entropy of the combustion products were computed at a com- 
bustion pressure of 300 pounds per square inch absolute (20.41 atm). 

Procedure for exit conditions . - Equilibrium composition, mean 
molecular weight, pressure, derivative of the logarithm of pressure with 
respect to the logarithm of density at constant entropy y , enthalpy of 

the products of combustion, specific heat at constant pressure, coeffi- 
cient of viscosity, and coefficient of thermal conductivity were computed 
for each equivalence ratio by assuming isentropic expansion for four 
assigned exit temperatures selected to cover the exit pressure range from 
the nozzle -throat pressure to about 0.1 atmosphere. 

Interpolation . - Parameters for pressures at and near the nozzle 
throat and for pressures corresponding to altitudes of 0, 10,000, 20,000, 
30,000, 40,000, and 50,000 feet were interpolated by means of cubic equa- 
tions between each pair of the assigned exit temperatures. The functions 
and their first derivatives used in the interpolations are described in 
reference 1. 
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The errors due to interpolation were checked for several cases. 

The values presented for all performance parameters appear to be correctly- 
interpolated or in error at most by one or two units in the last place 
tabulated. 

Formulas . - The formulas used in computing the various parameters 
are given in reference 1 and are summarized as follows: 

Specific impulse, lb-sec/lb 


I = 294.98 


h c - h e 


1000 


Throat area per unit flow rate, (sq ft)(sec)/lb 
(pressure in atm) 


S t /v = 


_ 1.3144 T t 


P t Mt a 


Characteristic velocity, ft/ sec 
(P c = 300 lb/sq in. abs) 


c* = g P c s t/ w 

= 1.3899X10 6 S t /w 


( 1 ) 


( 2 ) 


(3) 


Coefficient of thrust 


C F = Ig/c* = 32.174 i/c* 

Nozzle-exit area per -unit flow rate, 

(sq ft)(sec)/lb (pressure in atm) 


S e / W 


0.040853 T e 
P e A 1 . 


(4)- 


( 5 ) 


Ratio of nozzle-exit area to throat area 
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Specific heat at constant pressure, cal/(g)(°K) 


C P nMT 


E %«£>, + E “itHT', Y i - E 

L i i i 


(7) 


Derivative of logarithm of pressure with respect to logarithm of density 
at constant entropy 


E Pi »i 

Ts ‘ P (D A -1) 

Coefficient of viscosity, poise 


V p i 

L* (Hi/Mi) 

Coefficient of thermal conductivity, cal/(sec)(cm)(°K) 


( 8 ) 


(9) 


k = l ‘( c p + !l) (10) 

When composition is assumed to be frozen, equations (7) and (8) 
become 

Specific heat at constant pressure assuming frozen composition 
cal / (g)(°K) 


E n i (C °) 1 


( c ) — 

' P frozen nM 


(11) 


Derivative of logarithm of pressure with respect to logarithm of density 
at constant entropy assuming frozen composition 


(rJ 


( C ) 

' P 'frozen 


s' frozen ( } H \c v / frozen 

P frozen M 


( 12 ) 


The values of viscosity ancf thermal conductivity for mixtures of 
combustion gases calculated by means of equations (9) and (10 ) are only 
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approximate. When more reliable transport properties, for the various •. 
products of combustion become available, a more rigorous procedure for 
computing the properties of mixtures may also be justified. 


THEORETICAL PERFORMANCE DATA 

The calculated values of the various performance par ame ters for a 
combustion pressure of 300 pounds per square inch absolute and at exit 
pressures corresponding to altitudes of 0, 10,000, 20,000, 30,000, 

40,000, and 50,000 feet are given in table II for ten equivalence ratios. 
The values of pressure corresponding to the assigned altitudes were taken 
from reference 6. As an aid to engine design, the values of the param- 
eters within the rocket nozzle for 80, 90, 100, 110, and 120 percent of 
the throat pressure are tabulated in table III. Equilibrium composition, 
r_, specific heat at constant pressure, coefficient of viscosity, coeffi- 
cient of thermal conductivity, and mean molecular weight in the combus- 
tion chamber and at assigned exit temperatures are given in table IV-. 

The mole fraction of F 2 was always less than 0.00002 and therefore was 

not tabulated in table IV. 

Parameters . - The parameters are plotted in figures 1 to 9. Curves 
of specific impulse for the six altitudes are shown in figure 1 plotted 
against weight percent fuel. The maximum value of specific impulse for 
the sea-level curve is 311.5 pound- se c onds per pound at 24.1 percent of 
fuel by weight. 

The maximum values of specific impulse and the weight percentages 
at which they occur were obtained by numerical differentiation of the 
calculated values and are shown in figure 2 as functions of altitude. 

The maximum specific impulse increases 22 percent for a change in alti-' v 
tude from sea level to 50,000 feet. 

Curves of combustion-chamber temperature and nozzle-exit tempera- 
ture for the six altitudes are presented in figure 3 as functions of 
weight percent fuel. The maximum combustion temperature obtained was 
4310° K at 21.4 percent fuel by weight. The maximums of the exit tem- 
perature curves occur near the stoichiometric ratio. 

Characteristic velocity and coefficient of thrust are plotted in ’ 
figure 4 and ratios of the area at the nozzle exit to the area at the 
throat are shown in figure 5 as functions of weight percent fuel. 

Curves of mean 'molecular weight in. the combustion chamber and 
nozzle exit are plotted against weight percent fuel in figure 6. 
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Curves of specific heat at constant pressure, coefficient of vis- 
cosity, and coefficient of thermal conductivity for six pressures are 
plotted in figures 7 to 9 as functions of weight percent fuel. 

Frozen composition. - In order to compare data based on the assump- 
tions of equilibrium and frozen composition during the expansion process, 
several additional calculations were made assuming frozen composition. 
These are presented in the following table together with corresponding 
equilibrium data for the stoichiometric equivalence ratio and expansion 
to two altitudes: 


Parameters 

Altitude 

Sea level 

50,000 

feet 

Equili- 

brium 

Frozen 

Equili- 

brium 

Frozen 

I, lb -sec/lb 

311.0 

287.9 

379.2 

336.2 

c*, ft/sec 

7019 

6690 

7019 

6690 

Cp 

1.426 

1.385 

1.738 

1.617 

S e / S t 

3.908 

3.131 

18.71 

12.90 

T e , °K 

3113 

2026 

2130 

1122 

Me 

20.72 

19.10 

21.14 

19.10 


The percentage differences in these parameters for frozen and equilibrium 
composition are considerably higher for expansion to 50,000 feet than for 
expansion to sea level. 

For a combustion-chamber pressure of 300 pounds per square inch 
absolute and an exit pressure of 1 atmosphere, the values of maximum 
specific impulse are 311.5 pound-seconds per pound at 24.1 percent fuel 
by weight for equilibrium composition during expansion and 290.0 pound- 
seconds per pound at 25.7 percent fuel by weight for frozen composition 
during expansion. 

Chamber pressure effect . - According to NACA data for liquid hydra- 
zine with liquid fluorine, the parameters c*, Cp, and S e /S t are very 

nearly linear with the logarithm of chamber pressure for a fixed equiva- 
lence ratio and expansion ratio. For the stoichiometric equivalence 
ratio, increasing chamber pressure by a factor of 2 resulted in a change 
of +1.0 percent for c *, and changes of -0.1 percent for C-p and 

-1.0 percent for S e /St for an expansion ratio of 20.41$ and changes of 
-0.6 percent for Cp and -3.3 percent for S e /S-t for an expansion ratio 
of 326.6. It is expected that the values of c*, Cp, and S e /S^ given 
in this report for liquid ammonia with liquid fluorine for a chamber 
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pressure of 300 pounds per square inch absolute may be used at other 
chamber pressures with similar small differences. Greater precision can 
be obtained by additional performance computations for other chamber pres- 
sures. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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TABLE I. - PROPERTIES OF LIQUID PROPELLANTS 


Temperatures in superscripts, °C 







Propellant *■ 

' Properties 

Ammonia 

Fluorine 

Molecular weight, M 

17.032 

38.00 

Density, g/cc 

a 0.68“ 33,4 

b 1.54 _196 

Freezing point, °C 

c - 77.76 

c -217. 96 

Boiling point, °C 

c -33.43 

C -187. 92 

Viscosity, centipoises 

a 0.255" 33,5 


Enthalpy of formation at Boiling point 



from elements at 25 °C, AH^., kcal/mole 

d -17.14 

d -3.030 

Enthalpy of vaporization, AH, kcal/mole 

c 5.581“ 33 * 43 

c 1.51-187.92 

Enthalpy of fusion, AH, kcal/mole 

c 1.351 -77 * 76 

c 0.372" 217 * 96 


a Reference 3. 
^Reference 4. 
Reference 2. 
^Reference 5. 


TABLE II. - CALCULATED PERFORMANCE OF LIQUID AMMONIA WITH LIQUID FLUORINE 
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Based on F 2 density of 1.54 at -196° C and NH* density of 0.68 at -33. 
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Figure 1. - Theoretical specific impulse of liquid ammonia with liquid fluorine. Isentropic expan- 
sion, assuming equilibrium composition; combust ion- chamber pressure, 300 pounds per square inch 
absolute; exit pressure corresponding to altitude indicated. 
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Figure 2. -" Maximum theoretical specific impulse and corresponds 
ing weight percent of fuel in propellant of liquid ammonia with 
liquid fluorine. Isentropic expansion assuming equilibrium 
composition; combust ion- chamber pressure, 300 pounds per square 
inch absolute; exit pressure corresponding to altitude indicated 
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-Figure 4. - Theoretical characteristic velocity and coefficient of .thrust of liquid 
ammonia and liquid fluorine. Isentropic expansion assuming equilibrium composi- 
tion; combustion chamber pressure, 300 pounds per square inch absolute; exit 
pressure corresponding to altitude indicated. 
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Figure 5. - Theoretical ratios of nozzle-exit area to throat 
area of liquid ammonia with liquid fluorine. Isentropic 
expansion assuming equilibrium composition; combustion- 
chamber pressure, 300 pounds per square inch absolute; exit 
pressure corresponding to altitude indicated. 
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Fuel in propellant, percent by weight 

Figure 6. - Theoretical mean molecular weight in combustion chamber and at nozzle 
exit of liquid ammonia, with liquid fluorine. Isentropic expansion assuming 
equilibrium composition; combustion-chamber pressure, 300 pounds per square inch 
absolute; exit pressure corresponding to altitude indicated.: • 
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Figure 7. - Theoretical specific heat at constant pressure of combustion 
. products (including energy of dissociation) of liquid ammonia with liquid 
fluorine... Isentropic expansion assuming equilibrium compos it ion ; combustion 
chamber pressure,. 300 pounds per square inch absolute; exit pressures as 
indicated . 
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Figure 8. - Theoretical coefficient of viscosity of combustion products of liquid 
ammonia with liquid fluorine. Isentropic expansion assuming equilibrium composi- 
tion; combustion- chamber pressure, 300 pounds per square inch absolute; exit 
pressures as indicated. 
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Figure 9. - Theoretical coefficient of thermal conductivity of combustion 

products of liquid ammonia with liquid fluorine. Isentropic expansion assum- 
ing equilibrium composition; combust ion- chamber pressure, 300 pounds per 
square inch absolute; exit pressures as indicated. 
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